Quantifying Entanglement with Witness Operators 
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We present an unifying approach to the quantification of entanglement based on entanglement 
witnesses, which includes several already established entanglement measures such as the negativity, 
the concurrence and the robustness of entanglement. We then introduce an infmite family of new 
entanglement quantifiers, having as its limits the best separable approximation measure and the 
generalized robustness. Gaussian states, states with symmetry, states constrained to super-selection 
rules and states composed of indistinguishable partides are studied under the view of the witnessed 
entanglement. We derive new bounds to the fidelity of teleportation d m in, for the distillable en- 
tanglement Ed and for the entanglement of formation. A particular measure, the PPT-generalized 
robustness, stands out due to its easy calculability and provides sharper bounds to d m i n and Ed than 
the negativity in most of states. We illustrate our approach studying thermodynamical properties 
of entanglement in the Heisenberg XXX model. 

PACS numbers: 03.67.-a 
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I. INTRODUCTION 

In recent years entanglement has been recognized as 
a physical resource central to quantum information pro- 
cessing. As a result, a remar kable research effort has 
been devoted to classifying and quantifying it. The first 
achievement in this direction was the identification of the 
entropy of entanglement [lj , Ee , as the unique measure 
of entanglement for puré bipartite states in the asymp- 
totic limit. It was shown that that m copies of a pure 
state can be reversibly converted into n copies of \<f) 
by local operations and classical communication (LOCC) 
if, and only if, itiEe(\iP)) — nEE{\4>))- This reversibil- 
ity is lost however when one considers the more general 
picture of mixed states. In this case two different en- 
tanglement measures, associated with the formation and 
distillation processes, respectively, have to be taken into 
account. On one hand the entanglement cost, Ec{p) Q, 
is the minimal number of singlots necessary to create the 
state p by LOCC in the asymptotic regime. On the other, 
the distillable entanglement, Eo{p) is the maximum 
number of singlets that can be extracted by LOCC from 
p. Another important measure connected to asymptotic 
properties is the relative entropy of entanglement, Er 
It is related to how distinguishable an entangled state is 
from a separable one and gives bounds to Eq and Ep. 

The finite copy case is more complex and the entropic 
quantities considered above are not applicable anymore. 
For bipartite pure states, where the reversibility is al- 
ready lost, the minimum set of entanglement measures 
characterizing deterministic and probabilistic transfor- 
mations were derived j^QS- The mixed case, however, 
is known only for very restricted situations and remains 
mainly unsolved. 
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Another approach for the quantification of entangle- 
ment is to measure the usefulness of a state to perform 
a given quantum information task. For example, the 
maximal fidelity of teleportation achieved by single copy 
LOCC ;6j, the maximal secret-key rate attainable by lo- 
cal measurements in a cryptographic protocol and the 
capacity of dense coding despite not being equal to 
any of the measures discussed so far, are clearly the best 
quantifiers when one of these protocols is analyzed. 

Entanglement in multi-partite systems exhibits a much 
richer structure than the bipartite case and its study is 
even more challenging. Already in the pure three qubits 
case, there are two different manners for a state to be en- 
tangled, in the sense that there are states that cannot be 
converted, even with a certain probability, in each other 
From the measures considered above only Er is un- 
ambiguously defined to multi-partite systems, although 
it is not the only one. 

It is thus clear that entanglement is a highly complex 
phenomenon, which cannot be quantified by only one 
measure. Then, a natural way to measure it is to use 
any quantie which satisfies some particular properties, 
being the monotonicity under LOCC the most important 
0, Il0| . In this axiomàtic approach any measure which 
does not increase, on average, under LOCC, called an 
entanglement monotone |lfjj | , is a good measure of entan- 
glement and, conversely, any meaningful quantifier has 
to be an entanglement monotone, or at least has some 
sort of weaker monotonicity under LOCC. 

A closely related problem to the quantification problem 
is the characterization of entanglement. The very funda- 
mental question whether a given mixed state is entangled 
or not is extremely difficult, being actually NP-HARD 
[ïï| . A possible approach is then to consider sufficient cri- 
teria for entanglement, such as the Peres-Horodecki [ï^ 
and the alligment tests. Nonetheless, the strongest 
manner to characterize entanglement is using entangle- 
ment witnesses (EW) ^[ ^jj . They are Hermitian oper- 
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ators whose expectation value is positive in every sepa- 
rable state. Therefore, a negative expectation value in a 
measurement of a witness operator in an arbitrary state 
is a direct indication of entanglement in this state. Fur- 
thermore, it was shown by the Horodeckis that a state is 
entangled if, and only if, it is detected by an EW [Ï3 |. 
A great deal of research has been devoted to the study 
of EWs, yarying from the their classification and opti- 
mization |l5l Il6l to their use in the characterization 
of entanglement in important, even macroscopic, physi- 
cal systems 0, 0, |2(j . Also optimal set-ups for local 
measurements of witnesses |2ll |22| and experimental re- 
alizations of witnessing entanglement were realized |23|. 
In spite of the determination of EWs for all states being 
also computationally intractable [ÜJ, different methods 
from convex optimization theory can be applied to the 
problem, leading to efficient approximative procedures 
to determine and even optimize EWs for arbitrary states 

The main objective of this paper is to show that EWs 
can be very helpful also in the quantification of entan- 
glement. The first measure related to EWs was due to 
Bertlmann et al |27l | and was shown to be equal to the 
Hilbert-Schmidt distance from the set of separable states. 
Brandao and Vianna j2^| took another significant step in 
this direction, showing that a measure derived from op- 
timal EWs of the most studied group of witnesses so far, 
the group of EWs with unit trace, was in fact equal to 
the random robustness, which led to the establishment of 
properties still unknown for the later, such as its mono- 
tonicity under separable trace-preserving superoperators. 

Besides the obvious benefit of increasing the number 
of entanglement measures known, EWs based quanti- 
fiers are particular interesting due to the possibility of 
performing experimental measurements of them, which 
could be important to the extension of entanglement to 
other areas of physics, such as thermodynamics and sta- 
tistical mechanics. Moreover, despite being necessary in 
general a complete tomography of a state to the deter- 
mination of its degree of entanglement based on an EW 
measure, any EW provides a lower bound to it, even when 
no information at all about the state is available. 

The paper is structured as follows. In Sec. II we briefly 
review the bàsic properties of multi-partitc optimal en- 
tanglement witnesses. In Sec. III wc definc a class of 
entanglement measures based on EWs, which includes 
several important already known quantities such as the 
negativity and the concurrence, and introduce a new infi- 
nite family of entanglement monòtones having the gener- 
alized robustness and the best separable approximation 
measure as its limits. In Sec IV we present further prop- 
erties of the considered measures and relate them to the 
localizable entanglement. In Sec. V it is shown how 
the methods developed in the last years to the charac- 
terization of entanglement based on convex optimization 
can be used to calculate approximately a large number 
of measures based on EWs. In Sec. VI possible exten- 
sions of our approach to Gaussian states are discussed. 



In Sec. VII we consider how the measures and their cal- 
culation are modified in states with symmetries. In Sec. 
VIII the qüestions of the amount of entanglement and of 
nonlocality in the presence of a super-selection rule are 
answered from the perspective of the studied measures. 
In Sec. IX it is shown that the three most successful 
approaches to the quantification of entanglement in sys- 
tems of indistinguishable partides can be easily accessed 
from the EWs based quanties. In Sec. X and XI the 
qüestions of bounds on the teleportation distance and on 
the distillable entanglement of a given quantum state is 
review using our measures. It is shown that they provide 
sharper bounds than the negativity for the majority of 
states. In Sec. XII we derive lower bounds to the en- 
tanglement of formation with any EW. Possible applica- 
tions of the measures are exemplified in Sec. XIII, where 
the derivation of two thermodynamic equations of state 
which take into account entanglement is presented. Fi- 
nally, in Sec. XIV we summarize our results and discuss 
future perspectives. 

II. MULTIPARTITE SYSTEMS AND OPTIMAL 
ENTANGLEMENT WITNESSES 

We consider a system shared by N parties {Ai}f =l . 
Following [2^, we call a k— partite split a partition of 
the system into k < N sets {Si}^ =1 , where each may 
be composed of several original parties. Given a den- 
sity operator pi...k G B{H\ ® ... ® H^) (the Hilbert space 
of bounded operators acting on -Hi (8 ... ® iïfe) associ- 
ated with some k— partite split, we say that p±.,,k is a 
m-separable state if it is possible to find a convex de- 
composition for it such that in cach pure state term at 
most m parties are entangled among each other, but not 
with any member of the other group of n — m parties. 
For example, every 1-separable state, also called fully- 
separable, can be written as: 

Pi...k = y^Pi\ipi)i(j>i\ ® ■■■ ® \ipi)k{ti>i\- (1) 

i 

Another example is the 2-separable states of a 3-partite 
split given by: 

Pl:2:3=^T f PiPi, (2) 

i 

where each pi is separable with respect to at least one 
of the three possible partitions (A:BC, AB:C and AC:B). 
For each kind of separable state there is a different kind 
of entanglement associated to it. We will say that a state 
is (to + l)-partite entangled if it is not TO-separable. It 
is clear that if a state is m-separable it cannot be n- 
entangled for all n > m. 

It is possible to detect (m+l)-partite entanglement us- 
ing entanglement witnesses. In order to do that, consider 
the index set P = {1, 2, k}. Let P m be a subset of P 
which has at most to elements. Then W is a (to + 1)- 
partite entanglement witness if 
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, R m such that 



(3) 



uLi pr = p and prn^ m = {}- 



Equation (3) assures that the operator W is positivc 
for all m-separable states. Thus, as the subspace of in- 
separable density operators is convex and closed, a state 
p is (m + l)-entangled if and only if there is a Hermitian 
operator satisfying equation (3) such that Tr(W p) < 



Usually one is interest in a selected group of witnesses 
operators called optimal. Two different dcfinitions of op- 
timal entanglement witness (OEW) exist. The first, in- 
troduced by Lewenstein et al 0] , is based on how much 
cntangled states a given entanglement witness (EW) W 
is able to detect: W is optimal iff there is no other EW 
which detects all the states detected by W and some 
other states not detected by W . The second definition, 
due to Terhal 0] , establish the concept of OEW relative 
to a chosen entangled state p. The p— optimal entangle- 
ment witness W p is given by 

TríWpp) = min TriWp), (4) 
weM 

where Ai is the intersection of the set of entanglement 
witnesses, denoted by W, with some other set C such 
that Ai is compact 31] . Note that every p-OEW is also 
an OEW accordingly to the first definition, whereas the 
converse may not be true. 

A general expression for entanglement witnesses was 
presented in [3Ï|. Every EW acting on fc-partite Hilbert 
space can be written as: 



(5) 



where P and the QiS are positive semi-definite, e > 0, I 
is the identity operator and 7j is the partial transposition 
with respect to partie i. Note that even (m+lVpartite 
EWs can be written in the form of equation (5) [33]. An 
important class of EW is the decomposable entanglement 
witnesses (d-EW), which can always be written as: 



W = P 



(6) 



This class will be particularly important in our discus- 
sion, since the set of entangled states detected by d-EW 
is invariant under LOCC 13411. 



III. DEFINITIONS AND BASIC PROPERTIES 

In this section we show how p-optimal EWs can be 
used to quantify all the different kinds of multipartite en- 
tanglement. First, an unifying approach, which includes 



several important entanglement measures (EM), will be 
presented. Then we will consider a new infinite family of 
entanglement monòtones [ïc| . 

A general expression for the quantification of entangle- 
ment via EWs is defined as: 



E(p) = max{0, - min TrCWp)}, 



(7) 



where Ai = WC)C, and the set C is what distinguish the 
quantities. We call witnessed entanglement any measure 
expressed by equation (7). 

Some wcll known EM can be expressed as (7). The 
first, introduced by Bertlmann et al |27j . is: 

B(p)= max [min Tr (W a) -Tr(Wp)], (8) 

\\W— -T||2<1 cgS 

where W £ W. B{p) was shown to be monotonic decreas- 
ing under mixing enhancing maps |35| and to be equal to 
the 7í s -distance of p to the set S of fully-separable states: 



B(p) = D(p) 



mm 

o-6S 



(9) 



The second is the negativity, i.e., the sum of the neg- 
ative eigenvalues of p TA (the pa rtial transpose of p with 
respect to subsystem A) [3(| |33, H^] . It is easily seen that 
Aí can be written as: 



Aí(p) = max{0, 



min Tr(W TA p)\. 

o<w<i 



(10) 



Another quantie is the maximal fidelity of distillation 
under PPT-protocols, introduced by Rains |39j . 

Fd(p) = -,+ max{0, - min Tr(W TA p)}, (11) 
d weM 

where M = {W | (1 - d)I/d < W < I/d, < W Ta < 

21 /d} H3. 

The last is the celebrated Wootter's concurrence of two 
qubits, which can be written, accordingly to Verstraete 
El, as 



C(p) = max{0, 



min Tr((\A)(A\) T *p)}, (12) 

AeSL(2,C) 



where \A) denotes the unnormalized state (A(g)I)\I) with 
\I) = ^ I"), det(A) = 1. 

Assuming that the set C is also convex, which is the 
case of all the quantities considered in this paper, except 
the concurrence, it is possible to apply the concept of 
Lagrange duality from the theory of convex optimization 
to the problems represented by equation (7) 42]. Re- 
markably, the dual measures obtained are those related 
to mixin g p roperties, such as the robustness of entan- 
glement ÜJ], introduced by Vidal and Tarach, and the 
best separable approximation measure ^3], introduced 
by Karnas and Lewenstein. Moreover, since in all the 
cases considered here there always exist a strictly feasi- 
ble point, i.e, a W G relintAÍ (denoted in the convex 
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optimization literature by Slater condition), the optimal 
solution of the primal and dual problems are the same, 
i.e., the primal and dual measures are equal 

We now show that the dual representation of the gen- 
eralized robustness of entanglement Rg{p) HÜ, i-e., the 
minimal s such that 



P 



1 



(13) 



is separable, where ir is any, not necessarily separable, 
density matrix is given by (7) with M. = {W G W | W < 
I}. Following [42|. the Lagrangian of the problem is given 
by 



L(W,g(a),Z) = Tr(Wp) +Tr(WZ) 
-Tr{Z) - i g{a)Tr{Wa)da, 



(14) 



crES 



where Z, g(a) are the Lagrange multipliers associated 
with the constraints W < I and Tr(Wa) > Ver G S, 
respectively. Note that since the definition of EW is a 
composition of infinite constraints, its Lagrange multi- 
plier is a generalized function |4*Ej . The dual problem is 
then 

minimize Tr(Z) (15) 

subject to Z > 

g{a) > 0, Ver G S 

p + Z = J g(a)ada. 

Since g(a) > 0, the integral in the constraints above is 
a separable state. Conversely, any separable state a Q is 
obtained with the choice of g(<j) = 5{a — <j ). It is then 
easily seen that the result of (15) is the generalized ro- 
bustness. The dual representation of the robustness of 
entanglement, R(p), has, instead of W < I, the con- 
straint Tr(Wa) < 1, V a G S. 

The best separable approximation mcasure BSA(p) 
|44j is the minimum A such that there exist a separable 
state cr and a density operator Sp satisfying 



p = (1 - X)a + XSp. 



(16) 



It can been seen that the dual representation of BSA(p) 
is given by (7) with M = {W G W \ W > -I}. 

In |28j it was shown that the random robustness R r {p) 
431, i.e.. to the minimal s such that 



p + s(I/D) 



1 



(17) 



is separable, is equal to equation (7), with M = {W G 
W | Tr(W) = D}. This result can also be easily derived 
using the concept of Lagrange duality. 

In the next subsection we will introduce our new family 
of entanglement monòtones 



A. A New Family of Entanglement Monòtones 

If we let C be the set of Hermitian matrices W such 
that —ni < W < ml, where n, m > 0, then the quantity 
derived from (7) will be denoted by E n:m . 

Proposition 1. E n:m is an entanglement monotone for 
every n, m > 0, i.e. 



PiEn-.jnip'i) < E n:m (p) 



(18) 



where p\ is the final state conditional on the oceurrence 
of the classical variable i, which oceurs with probability 
Pi at the end of a LOCC protocol. 

Proof. It suffices to consider final states of the form 

P \ = AtpAl/pi (19) 

with pi = Tr[AipA\], where the Kraus operators 
Ai, A M are given by Ai = A\ ® ... ® A\ and satisfy 

Y, i PiEw(pd = EiPi^{0,-Tr(W p ,p' i )} 
- J2 k -Tr{A\W p , h A k p) < -Tr{W pP ) = E w (p) [ Uj 

where k sums only the terms such that 
max{0, — Tr(W p 'p' i )} is different from zero. In the 

last inequality we used that W = J^k ^kS^p' ^-k < 
mY, h A{A k < ml, W = £ fc 4woA* > 



•'E k A{A k > 



-ni, and that W„ is optimal. 



□ 



Note that the proof of proposition (2), with minors 
modifications, applies also to Aí and F4. 
The dual representation of E m:n (p) is 

minimize ms + nt (21) 

subject to p + S7ri = (1 + s — t)a + Í7T2 

where ttí are density matrices, a is a separable state and 
s, t > 0. From (24) we find that 

lim E n:m (p) = nBSA(p), lim E n:m (p) = mR G (p) 

m — >oo n — >oo 

Actually, the equalities above are already vàlid when one 
of the numbers is sufficiently larger than the other. The 
elements of this new family of EMs can be interpreted as 
intermediate measures between the generalized robust- 
ness and the best separable approximation. Note that 
for every distinct rational number n/m within a cer- 
tain finite interval, the E m:n are genuine different EMs, 
meaning that there is no positive number c such that 
E m:n = cE m i. m , ïínjm^ rí /m'. 

If we consider that C is the intersection of set of Hermi- 
tian matrices W such that —ni < W < ml with the set 
of decomposable entanglement witnesses, a new family of 
entanglement monòtones, denoted by E^.^ is defined. 
To see that they are indeed EMs, all we have to note is 
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that for every A i = A\®...®A k i , A\WAi is a decompos- 
able EW whenever W is. Therefore, proposition (2) also 
applies to them. 

It is possible to derive several other familics of EMs 
considering intersections of the sets C of diffcrent entan- 
glement measures which can be written as equation (7), 
such as those given by equations (10-12). 



IV. MULTIPARTITE ENTANGLEMENT 
HIERARCHY 

We now discuss more about the different kinds of mul- 
tipartite entanglement introduced in the second section. 
Usually the set of separable states is regarded to be com- 
posed of all sates which can be created by LOCC proto- 
cols. In this sense, given a specific split and considering 
that each part of the split can perform global quantum 
operations on its subsystems, only 1-separable states can 
be properly identified as separable. However, one might 
also be interest in the situation where some of the par- 
ties are allowed to perform join operations. In this case, 
the different types of entanglement play an important 
role. Consider, for example, the situation where k par- 
ties want to create a common quantum state and each 
one is connected to the others via a quantum channel. 
If they all agree in using their channels, every state can 
be prepared and the situation becomes trivial. However, 
suppose that they agree that only m < k parties will use 
their quantum channels, where the probabilities of which 
parties will be involved are given by pi. At the end of the 
protocol they will share an ensemble of states, {pi,Pi}, 
which clearly does not have m + 1-partite entanglemnt. 
Now, since erasing classical information cannot create en- 
tanglement, we are lead to consider the different kinds of 
entanglement discussed before. This property is reflected 
in the condition that every goog entanglement measure 
should be convex, which we show for every quantity de- 
fined according to equation (7) 

Proposition 2. E is a convex function for any choice 
ofC, i.e., 



E £ 



pí pí 



< 



(22) 



whenever the pi are Hermitian, and pi > with Y jí Pí = 



1. 



Proposition (3) follows from the convexity and the con- 
cavity of the max and min functions, respectively. 

Consider a given fc-partite split of a multi-partite sys- 
tem p. It is possible to attribute (k — 1) numbers, E m , 
1 < m < k — 1, where each one quantifies one type of 
multi-partite entanglement of the system. It is easy to 
see from equation (3) that all constraints imposed to an 
EW which deteets m-partite entanglement (m-EW), are 
also imposed to every n-EW, with n > m. Hence, the 



following order between the E m holds: 

E m {p) >E n {p), V n > m. 



(23) 



E 1 (p), formed by the OEW with respect to the fully 
separable states is an upper-bound to all other E(p), 
including those with respect to other splits formed by 
grouping several original parties into one. This means, 
for example, that in a 3-split, E 1 (p) is greater or equal 
to the bipartite entanglement of any of the three 2-splits, 
namely A-BC, AB-C and AC-B. Actually, it is possible to 
establish a complete hierarchy in the proposed measures 

An interesting measure of entanglement for multi- 
partite systems is the localizable entanglement, intro- 
duced by Verstraete et al |50j| . Given a quantum system 
of n parties p, the localizable entanglement Eij (p) is the 
maximal amount of entanglement that can be created, 
on average, between the parties i and j by performing 
a single-copy LOCC protocol in the system |5l). More 
specifically, if at the end of a LOCC protocol we have an 
ensemble of states p = {pk , p\ 3 } , where pu is the proba- 
bility that the reduced state of the parties i and j is p % i , 
the LE is then given by 



E VJ 



max 



^2PkE(p 



(24) 



where E{p) represents, in this paper, one measure based 
on OEWs. The LE has the operational meaning which 
applies to situations in which out of some multipartite 
entangled state one would like to concentrate as much 
entanglement as possible in two particular parties |50| . 
which could be used later, for instance, in some quantum 
information task. 

Proposition 3. Consider a multi-partite state p. Then 
E& m <K m (p)> V í,i,n,m. (25) 

Proof. As in the proof of proposition 2, it sufhces to con- 
sider final states of the form 



pf =T r/ij {A w A\/ Pl ) 



(26) 



withp; = Tr[AipAj], where Tr/y stands for the partial 
trace of all parties, except i and j. The Kraus operators 
Ai,..., Am are given by A\ = A] ® ... ® A\ and satisfy 

*%,m = EiPiEnUpí 3 ) = Ei max{0, -Tr(I ® W# Pl )} 
- E fc -Tr{A\l <g> WuA kP ) < -Tr(W pP ) = E n:m (p) 



Pk 

where fc sums only the terms such that 
max{0, — Tr(W p ' i p' l )} is different from zero. 
In the last inequality we used that the EW 
W = Y^Ail ® WÜA* : m 



W = J2 k A\l ® W»A k > -nJ2 k AlA k > -ni and that 
W p is optimal. Note that proposition (5) also applies to 



E, 



PPT 



□ 
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The following relation between the negativity, Af(p), 

Proposition 4. 

M{p)<E^{p)<dN{p) (28) 
Proof. For every positive operator M, we have 

Xmax (M Ta ) < Kaax (M) < d\ max (M Ta ) (29) 

where the first (second) inequality is saturated for sep- 
arable (singlet) states. Hence, as < W < 1 implics 
W Ta < 1, we find 



N(p) = - min Tr(W TA p) 
o<w<i 

<- min Tr(W TA p) = E^{p) 

W T A<I 

W>0 



(30) 



where we have used that the optimal decomposable EW 
for a bipartite system has always the form W Ta , W > 0. 
From equation (32) we also find that W > 0, W Ta < I 
implies < W < dl. Thus, 



E^ P ?{p) = - min Tr(W TA p) 
w Ta <i 

W>0 

< - min TríW TA p) = dAÍ(p) 

0<W<dI 



(31) 



□ 



The second inequality is strict for example on the state 



(32) 



I-d(P+) T * 
P= <P-d 

where P + is the maximal d x d entangled state. 

V. NUMERICAL CALCULATION 



The lack of an operational procedure to calculate en- 
tanglement measures in general is ultimately related to 
the complexity of distinguish entangled from separable 
mixed states, which was shown to be NP-HARD [ïlj. 
Since an operational measure, which has positive value in 
every entangled state, would also be a necessary and suf- 
ficient test for separability, we should not expect to find 
one. Nonetheless, some approximative numerical meth- 
ods based on convex optimization have been proposed to 
the separability problem |2J, |2jJ |26( • What we will show 
in this section is that these methods can also be used to 
calculate, approximately, the witnessed entanglement. 

The first one, proposed in |2 1| by Brandao and Vianna, 
linked the optimization of EWs with a class of convex 
optimization problems known as robust semidefinite pro- 
grams (RSDP). Although RSDPs belong to NP-HARD, 
some well known probabilistic relaxations, which trans- 
forms the problem in a semidefinite program (SDP), 



werc applied, leading to a method of optimizing pseudo- 
EWs (operators which are positive in almost all separable 
states) to every multipartite state and with respect to all 
types of entanglement. 

The second approach, due to Doherty et al [25l l34| , 
was actually the first method to the separability prob- 
lem based on SDP. Using the existence of symmetric ex- 
tensions for separable states and the concept of duality 
in convex optimization, a hierarchy of SDPs, where the 
(k + \) th test is at least as powerful as (k) th (but de- 
mands more computational cffort), was constructed to 
detect entanglement. In each step k, an OEW with re- 
spect to a restricted set of EWs, which converges to the 
whole set of EWs in the limit of k — > oo, is obtained. This 
method can be used, however, only for the entanglement 
with respect to the fully-separable states, E 1 . Note that 
the further constraints that we demand to the EWs can 
be incorporated in the SDP, since they are linear matrix 
(in)equalities. 

The last method, introduced by Eisert et al [2(|, is 
based on recently developed relaxations of non-convex 
polynomial problems of degree thrcc in a hierarchy of 
SDPs, which converges to the solution of the original 
problem as the dimension of the SDP reaches the infinity 
|52|. One of the applications of this method is the min- 
imization of the expectation vàlues of EWs with respect 
to pure product states. Therefore, it can be used together 
with the second method discussed to lower the value of 
Tr(Wp), where W is a non-optimal EW determined by 
some step of the hierarchy. 

We would like to stress the complementary charac- 
ter of these methods. Whereas the first method usu- 
ally provides upper bounds to E{p) 1 since it only deter- 
mines pseudo-EWS, the second and third provides lower 
bounds to E(p), as the EWs resulting from them are 
non-optimal. Although only the first one can calculate 
E m (p), for m > 1, the number of constrains imposed 
grows exponentially with m. Thus, in most of cases, we 
will restrict ourselves to the determination of E , which 
is an upper bound to all other types of multi-partite en- 
tanglement (see section IV). 

Note that measures restricted to decomposable EWs 
can always be exactly calculated, in the worst case, by a 
semidefinite program. 



A. Example I 

As a first example we calculated E\. x and É* :1 , < 
n < 1, for the following family of states 

p q = q\W){W\ + {l-q)\GHZ){GHZ\, Q<q<\ (33) 

where \W) = (|001) + |010) + |100))/V3 and \GHZ) = 
(|000) + |lll>)/v / 2- The results are plotted in figures (1) 
and (2). When n « 1, £^ :1 = nBSA and we see that for 
all q there is not product vectors and even biseparable 
vectors in the range of p q . In the other limit, where 
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|[)00>+|0]l>+|l02>l·|ll3> 

(1) 



jooo>+ j^t|o n >+| i o i >)+j 1 1 : 

(2) 



|Í·)0()>+|()II>+|1I2> 



|ooo>+[m> 



|(l(M>+fí) [()>+[ 100> 

(6) 



FIG. 3: Representative states of the five distinct classes of 
3-entangled states. 



, 



FIG. 1: (Color online) El :1 (p q ) for < n < 4 and < q < 1. 





05 1 1.5 2 2.5 3 3.5 



FIG. 4: E^.^pq), < n < 4.5, for the states (1-5) of figure 
(3). 



FIG. 2: (Color online) E% :1 (p q ) for < n < 4 and < q < 1. 



E} rl — Rg, we find that the generalized robustness of 
cntanglement with respect to biseparable states is the 
same for all p q with q < 0.7. 



C. Example III 



As a final example, we present a numerical compari- 
son between Af and E^j = Rg PT ■ The bipartite PPT- 
generalized robustness can be determined as easily as the 
negativity. Actually, it can be written as 



B. Example II 



R P G PT {P) = 



(34) 



The classes of entangled states equivalent by SLOCC 
for 2 X 2 X n systems were determined in [53( and can 
be represented by the states (1-5) of figure (3). The 
arrows indicate which transformations are probabilistic 
possibles. E^a , < n < 1, was calculated for each 
of these and plotted in figure (4). Note that for all n 
considered, the incomparable states (2-3) and (4-5) with 
respect to state transformations have approximately the 
same E^^. 



where X m ax (P Ta ) is the maximum eigenvalue of the par- 
tial transposed projector onto the negative eigenspace of 
p T A We have generated 10 5 random states using the al- 
gorithm presented in |38| and plotted in figures (5) and 
(6). 

Although dj\í > Rq PT > Aí (see section IV), we see 
from figures 5 and 6 that Rq PT < 2AA for the majority 
of states. 



8 



displacement [54j. The symplectic matrix is given by 




FIG. 5: (Color online) Ra PT x N for 10 5 4 x 4 random sates. 




FIG. 6: (Color online) Rc PT x N for 10 5 6 x 6 random sates. 



VI. GAUSSIAN STATES 

We considered n distinguishable infinite dimensional 
subsystems, each with local Hilbert space H — £ 2 (TZ n ). 
A Gaussian state is characterized by a density operator 
whose characteristic function Xp( x ) = Tr[/9W(a:)] is a 
Gaussian function |54j . We can write, for every Gaussian 
state p, 

p = n- n í dxe-* x ''i x+idTx W(x), (35) 
Jn 2n 

where W(x) = exp[— íx T R] are the displacement opera- 
tors and R = (X u P u X 2 , P n ), with [X k ,Pi] = i5 kl . 
The matrix 7 > i J n is a 2n x 2n real matrix called cor- 
relation matrix (CM) and d is an 2n real vector called 



J„ = 0Ji, Ji = 



k=l 



-1 

1 



(36) 



Note that the displacement of a state can always be ad- 
justed to d = by a sequence of unitàries applied to 
individual modes. This implies that d is irrelevant for 
the study of entanglement. Thus, we set d = for now 
on without loss of generality. 

The optimization of EWs for states of infinite dimen- 
sion is completely infeasible. Nonetheless, we can still 
obtain meaningful quanties if we restrict it to a simpler, 
but sufficiently large, set of operators. An obvious choice 
would be the restriction to Gaussian entanglement wit- 
nesses (GEW), i.e., Gaussian operators which are positive 
in separable Gaussian states. Unfortunately, none Gaus- 
sian entangled state is detected by a GEW. Assuming 
that G is a GEW with covariance matrix T, we find 



Tr(pG) = í 



dxe~i xT ^ x+c > 0. 



(37) 



Anothcr possible class of operators then is given by 

Wç = {Q EB{H®H) \Q = 2 n I-G}, (38) 

where G is a Gaussian operator and / the identity oper- 
ator |5ï|. In the next proposition we show that E^. m (p) 
given by 

Eg.^p) =max{0,- min Tr(Qp)} (39) 
Qew e , Q<I 

can be very efficiently numerically calculated by a simple 
semidefinite program. 

Proposition 5. 

Eg :m (p) = max{0, det(T + 7 )" 1/2 - 2"} (40) 

where the matrix T € M2 n {TV) is obtained by the following 
SDP determinant maximization problem 

max dct(r + 7)~ 1/2 



subject to 



1 < r < 1, 



1 12 



r > 



r 2 + ítJ ± 12 

fi2 fi-jj 



> 



J n + r A \ 
A D(A) ) 



>o, 



Sfc-1,2/-1 Sk,l 
Sk,l Sfe-1,2/ 



(41) 



> 



flfe-1,21-1 > 0, a fc _i,aí > 0, s kl > 2, k = l.J, i = 1...2 1 '" 

where A is a n x n lower triangular matrix comprised 
of additional variables, D(A) is a diagonal matrix with 
same diagonal entries as those of A, l is the smallest 
number such that 2 l > n, and sqj = A^ if 1 < i < n and 
sq i — 2 if n < i < 2 l . 
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Proof. Consider the following structure for the bipartite 
Gaussian operator G 



G= / dxe-i* Tx W(x), 



(42) 



where T T — T > £ M2 n (7?.), with modes 1 to m and 
rn+ 1 to n belonging to Alice and Bob, respectively. The 
optimization objective Tr{Qp) 1 where Q — 2 n I — G, can 
be written as 



2 n -7T-" / dxe-* xT{T+l)x = 2"-dct(r+ 7 )- 



1/2 



(43) 



From the Jamiolkowski isomorphism, Q is an EW iff 
the map Q defined as Q = J® Q(P+) = 2"/ - 1® £(P + ) 
|5q is positive, which is equivalent to p' = I®G{p) < 2™/, 
for every density operator p. The covariance matrix of 
p', 7', can be written as |S3 

7' = S^VS, (44) 

where S € Sp(2n, 1Z) and <r is the covariance matrix 

o- = diag(/ii,/ii, (in, fin) 

corresponding to a tensor product of states diagonal in 
the number basis given by 



M 1 



fc=0 



1 



\k)u(k\, (45) 



| fc)i being the k-th. number state of the Fock space Hi 
[58j . The symplectic transformation (47) is reflected in 
the Hilbert space level by an unitary transformation: 
G = U G' U (S) . Since we are considering bounded 
operators, the pi must be non-negative. We also see that 
the positiveness of Q is equivalent to 



(46) 



Thus, since we are only considering Gaussian operators, 
equation (49) is satisfied iff 7' > for every 7 > iJ, 
where 7' and 7 are the covariance matrices of p' and p, 
respectively. Following Giedke and Cirac [Ü]], one finds 
that Q is positive iff 

min max {z^{M + iJ)z,z^(M-iJ)z}>0 (47) 

where M = f 2 -f? 2 (fi)- 1 f 12 and f = (I® A)T(I® A), 
with A = diag(l, —1, 1 , —1, —1). The matrices Vi are 
such that 



Ti r 12 

2 



r i2 r 2 



We now express condition (50) as a linear matrix in- 
cquality. Equation (50) is equivalent to 

z\M + iJ)z>Q Vz6C 2 " s.t. z\M -U)z <0 
z f (M-iJ)z>0 VzeC 2 " s.t. z f (Af + iJ)z<0 



An important theorem of matrix analysis, known as S- 
procedure, can be stated as follows: a quadratic function 
in the variable x, G(x), is positive for all x such that 
H{x) > 0, where H{x) is another quadratic function, iff 
there exists a positive real number t such that G — tH > 
0, for all x 59]. 

Applying it to the two conditions above we find that 
equation (50) holds iff there exists a positive number r 
such that 



M + tíJ>0, -Kr<l 



(48) 



We now use another fact of matrix analysis which says 
that the constraints on the Schur complement R > 
0, Q - SR^S 11 > and ker(i?) Ç ker(Q) are equiv- 
alent to 



Q 
s T 



> 



Hence, applying it to equation (51), we find that a Gaus- 
sian operator Q is an EW iff there exists a real number 
— f < t < 1 such that equation (44) holds. 

From the Williason decomposition, we see that Q < I 
is equivalent to 



n 



Vi 



= 2" det(7 + (S T )- l TS- 1 )- í < l , (49) 



where r\i ate the symplectic eigenvalues of Q. Since S 
is symplectic, one has S T J n S — J„, so that det(7 + 
(S^TS-^ = det[jJ(^)- 1 (J„ + r)5- 1 ]- 1 = 
det(J„ + r) 1 . The proposition then follows from ref- 
erence 60], which presents a LMI representation for the 
inequality det(A) 1 /" 1 > t, where A is a positive m x m 
real matrix. □ 



VII. STATES WITH SYMMETRY 

Entanglement measures have usually their calculation 
greatly simplified when the state in question has certain 
symmetries. Following [6f|], let G be a closed group of 
product unitary operators of the form U — U\ ® Ui. 
Dcfining the projection 



P(A) = j dU UAU*, 



(50) 



for any operator A on H\ ® H2 |63j , where dU is the Haar 
measure of G, we say that an operator M is invariant 
under G if P(M) = M, which is equivalent to [U,M] 
( i for all U E G. Consider now the determination of 
any measure expressed by equation (7). If the state in 
question has the property P(p) = p, then one can restrict 
the optimization in (7) to operators with the this same 
symmetry. More specifically, 



E(p)=- min Tr(Wp) 
= - min TrCP(W)p) 



min Tr(WV(p)) 

weM 

min Tr(Wp), (51) 
weP(M) 
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where P(M) = {W e M \ W = P(W)} |S3 . 

As an example, consider the isotropic states p p on C ( 

C d 



p p =pP+ + (l-p) 



d 2 ' 



(52) 



where P + = |<E> + )(<I> + | is the maximally entangled state. 
It can be shown that P + and the identity are the only 
operators which commute with all unitàries of the form 
U ®U* . Hence, the OEWs for p p can be written as 



W( Pp ) = a{p)P+ + b(p)I. 



(53) 



Since {ip\P + \^) < í/d, for every separable state \ip), we 
find 

P t s í(n + l) P + (1 -^" +1) -l n<d-L.,. 



l 



n > d. 



As the OEWs for this family of states are decomposable, 
equation (57) is also vàlid to E^.f T . 



VIII. SUPERSELECTION RULES 

The effect of superselection rules (SSR) in theory of 
cntanglement has been studied recently under a number 
of different perspectives [6J, |63, lMJ l67| . Two striking 
features emerge from the existence of a SSR. The entan- 
glem ent of a given state under SSR is usually reduced 
[65j and the notion of nonlocality has to be redefined, as 
there exists separable states that cannot be created by 
LOCC [Ü). In this section we show how the witnessed 
entanglemet fit in each of these scenarios. 

Following Bartlett and Wiseman j6^, we define a SSR 
as a restriction on the allowed local operations on a sys- 
tem, associated with a group of physical transformations 
G. An operation O is G-covariant if 



0[T{g)pT\g)] = T (g)0[p]T\g) , 



(55) 



for all group elements g S G and all density operators p. 
Then the SSR associated to G is the restriction on the al- 
lowed operations on the system to those G-invariant. As 
these restrictions make a state p indistinguishable from 
the states T{g)pT^(p) for all g £ G, it is convenient to 
describe p by the G-invariant state 



G(p) = / dgT(g)pTHg), 



(56) 



they can produce from G(p) by unconstrained LOCC. If 
E is an entanglement monotone, any LOCC applied to 
G(p), can, on average, at most maintain E(G[p])- Since 
it is always possible to reach this bound applying local 
swap operators, we have that the maximal amount of en- 
tanglement produced by SSR is exactly E(G[p]). Hence, 
from section (VII), it follows that, under a G-SSR 



E(p) = max{0, 



min Tr(W p)\, 

W£Ç[M] 



(57) 



where Q[M] = {W e M \ W = G[W}} [6g. 

We now consider the effect of SSRs in the notion of 
locality. The states that can be prepared locally in the 
presence of a G-SSR are those which can be written as 
(1), with each \ipi)k being G-invariant. It is possible 
to detect nonlocal states with witness operators, defin- 
ing a G-nonlocality witness (GW) as a Hermitian op- 
erator which satisfies equation (3), with |·0)pmpm (ip\ = 
G\\ip)p m p m (tp\], for all i and k. This nonlocal character 
of some states in the presence of a SSR can be quantified 
|67j . We can then, as it was done with entanglement, use 
GWs to perform this quantification. A witnessed nonlo- 
cality measure, Ng, will be any quantie given by equation 
(7), with the set of EWs substituted by the set of non- 
locality witnesses. It is easy to see that all properties 
discussed for E are vàlid for Ng- 

As an example, considered the following state 



p = i(|0) A <0|®|0) B <0| + |l)x<l|®|l)fl{l|) 



1^ 



I AB 



(58) 



where \V + ) A b = (\0)a\ p3 + M a\1) b) / VI 2 ) ■ Verstraete 
and Cirac have shown [64J that, although this state has 
a separable decomposition, it is not local when a particle 
number SSR exists, since all possible separable decom- 
positions have local states involving superpositions of a 
different number of partides. Any Hermitian matrix with 
positive diagonal entries is a G-nonlocality witness in this 
case. This should be contrasted with the case of a general 
EW, where an infinite number of inequalitics arc ncces- 
sary for its characterization. Calculating, for example, 
the nonlocal measure equivalent to i?oo:i = Rg, 

N G {p) = max{0,- min Tr(Gp)}, 

where Q = {G \ G tl > 0, G < /}, we find N G )(p) = 1/2, 
with G = -|01)(10| - |10)(01|. 



where dg is the group-invariant Haar measure. For mul- 
tipartite systems, where the SSRs are local, we have 
G[p] = G <S> ■■■ (8 G[p\- As it was shown in 65], the maxi- 
mal amount of entanglement which can be produced by 
LOCC in a register shared by all the parties, initially in 
a product state and not subjected to SSRs, from a state 
p, constrained by a G-SSR, is given by the entanglement 



IX. INDISTINGUISHABLE PARTICLES 

The study of entanglement in systems of indistinguish- 
able partides has been the subject of recent controversy 
© SEE El El- At least three different approaches to 
the problems have been proposed, namely, the entangle- 
ment of modes |69|| , the quantum correlations |72j | and the 
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entanglement of partides |73j . Each of these has its own 
advantages and drawbacks, and no consensus has been 
reached on which one is the most suitable. In this sec- 
tion we show how the proposed measures based on EW 
can be used to quantify entanglement in each of the three 
methods. 

We start with the entanglement of modes, proposed 
by Zanardi 69] , which suggests that the entanglement of 
indistinguishable partides should be calculated by any 
regular entanglement measure, using the density matrix 
in the mode-occupation, or Fock, representation. In this 
case it is clear that the determination of the witnessed 
entanglement follows straightforwardly. 

The quantum correlations, introduced by Schlicmann 
et al |Z2j, is motivated by the believe that no quantum 
correlations due to symmetrization (for bosons) or anti- 
symmetrization (for fermions) should be considered as 
true entanglement. Then, the characterization of entan- 
glement, for pure states, is determined by the Slater rank 
of the state, as opposed to the Schimdt rank usually con- 
sidered in distinguishable partides. Furthermore, Schlic- 
mann et al have shown that the concept of entanglement 
witness is also applicable to multipartite systems of in- 
distinguishable partides [ïil and, thus, the witnessed en- 
tanglement are well defined in this case too. 

The last approach, due to Wiseman and Vaccaro [73[ . 
is probably the best motivated one. The entanglement 
of partides is defined as the maximal amount of entan- 
glement, computed by a Standard measure, which Alice 
and Bob can produce between a quantum register, shared 
by them, composed of distinguishable partides by local 
operations. The amount of entanglement will clearly de- 
pends on the physical constraints imposed, which are in 
most of cases expressed as a SSR. Therefore the approach 
presented in the previous section to SSR can also be ap- 
plied in this case. 



X. TELEPORTATION DISTANCE 

In this section we derive lower bounds to the teleporta- 
tion distance, using E n:rn . We consider a quantum state 
p shared by k parties and ask what is the best possi- 
ble teleportation distance attained by a LOCC protocol 
when the parties form two groups and teleport a quan- 
tum state from one group to the other. 

Consider a teleportation protocol where a bipartite 
state pab is used as a quantum channel between Alice 
and Bob. Following the approach of Vidal and Werner 
to the negativity |3jj, we will first consider the single 
distance of a bipartite state defined as: 



A(P + ,p) = mï\\P + -P(p)\\ 1 , 



(59) 



where P+ is the maximally entangled state and the infi- 
mum is taken over LOCC protocols P. Using the convex- 
ity and the invariance under unitary transformations in 
the two terms of the absolute distance, and the invariance 



of P + under unitary transformations of the form U ®U* , 
we may assume that the optimal state which minimizcs 
equation (62), P opt {p), has undergone a twirling opera- 
tion |75j and, therefore, is a noise singlet, 



p p =pP+ + (1 -p) : 



f/2 



(60) 



The absolute distance of p p is given by ||P+ — p p \\i 
2(1 - p)(d 2 - l)/d 2 . From equation (57), 



|P+-p„||i = 2(l 



rpPPT 



(Pp) 



)■ 



(61) 



From the monotonicity of E^-i T under LOCC, we find 
that, for n > d, 



Proposition 6. 

A(P+,p)>2(l 



(P) 



)• 



(62) 



Since E^.\ T (p p ) — 2AÍ(p p ), for n > d, we see that 
EnK T provides, when E^ T < 2Af, a sharper bound 
than the one derived from the negativity. In the limit 
case n — > oo already, where E^ T is equal to the PPT- 
generalized robustness, we see from section (V) that the 
new bound is indeed sharper for the majority of states. 

A measure of the degree of performance of a quantum 
channel is the teleportation distance 



d(A)= / #||0-A(^)||i 



(63) 



As it was shown by the Horedecki family Q , the minimal 
teleportation distance that can be achieved when using 
the bipartite state p to construct an arbitrary teleporta- 
tion channel is given by 

d mm (p) = -^—A{P+,p). (64) 



1 



Therefore, 



dminip) < -7^7(1 ~ i±^§™), n > d. (65) 
d + 1 d 

Until now we have just adapted Vidal and Werner's rea- 
soning for the negativity to the E^.f T . Nevertheless, as 
opposed to N , E^.f T are also defined to multi-partite 
systems. 

Proposition 7. Consider a quantum state p shared by 
k parties. Let p 1 -- m -( m + 1 )-- k denote a bipartite split of the 
system, where the parties 1 to m and m + 1 to k form 
two groups. Then, 

j 1 A..m:{m+Ï)..k\ ^-D ^ 1 + (j^T T ) 1 (p) ^ 
" ,n.m\p J -i n , ,U jj )i 



D + l 



(66) 

V 1 < m < k, where D stands for the minimum of the 
dimensions of the two groups. 

Proposition (8) follows from the upper bound to all 
types of entanglement provided by E . Equation (57) 
is saturated, for example, on the fc-partite GHZ state 
\Vghz) = 1A/2(|00...0) + 
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XI. UPPER BOUNDS FOR THE DISTILLABLE 
ENTANGLEMENT 

We now move on to show another application of the 
family E± :m , namely bounds to the distillable entangle- 
ment of bipartite mixed states. We first derive the fol- 
lowing additivity property 

Proposition 8. 

E n a{p®p) <E n:1 {p) 2 + 2E n .. 1 {p), Vn>l. (67) 

Proof. Consider the dual representation (24) of E n:í . Let 
s, t, cr, tti and 7T2 be variables which minimize (24). Then 
we find E n: i(p) — s + nt, with 

p = (1 + s - t)<7 + ítt 2 - stt 2 . (68) 

Thus, 

p<gp=(l + s- t) 2 a <g a + í(l + s - t)a <g 7r 2 
-s(l + s - t)a <g 7T1 + í(l + s - í)tt 2 ® CT + Í 2 7T 2 ® 7T 2 

— SÍ7T 2 ® 7Tl — S(l + ,S — í)"7Tl ® (7 — SÍ7Tl ® 7T 2 + S 2 7Tl ® 7Tl 

= [1 + (2s + s 2 + í 2 ) - (2t + 2st)]a <g> cr 

+ [í(/9 (8) 7T 2 + 7T 2 (g) p) + SÍ(7Tl <g) 7T 2 + 7T 2 (g) 7Tl)] 
-[s(p (g) 7Tl + 7Tl (g) p) + S 2 7Tl (g 7T1 + Í 2 7T 2 (g 7T 2 ], 

where in the last two lines we used that 
1 



1 + s-t 



(p + STTl - Í7T 2 ) . 



It is therefore easily seen that if E n: i(p <g p) = s' + ni! 1 
then s' + nt' < s 2 +t 2 + 2s + ra(2í + 2sí). Hence, as n > 1, 

£;„ : l(p) 2 + 2-B„ : i(p) - ^ n: i(/>® p) 

= s 2 + n 2 t 2 + 2nst + 2s + 2nt - s' - ni! 
> s 2 + n 2 i 2 + 2nst + 2s + 2nt - s 2 - t 2 - 2s - 2nt - 2nst 

= i 2 (n 2 - 1) > 0. 

□ 

We can define a family of quantities close related to 
E n -.i by 



LE n:1 (p)=log 2 (l + E n:1 (p)). 



(69) 



The LE n: i(p) are non-increasing under trace preserving 
separable operations. From proposition (9) we find that 
they are also weakly-subadditive. Indeed, for n > 1, 

LE n:1 (p®p) <\og 2 ((l + E n:1 (p)) 2 ) = 2LE n:1 (p). (70) 

Note that the same results are also vàlid to E^ T . We 
now can state the main result of this section 



Proposition 9. 

E D (p) < LE n:1 {p), Vn> 1. 



(71) 



where Eu(p) is the distillable entanglement of the bipar- 
tite state p. 



Proof. The proof of proposition (10) is basically an appli- 
cation of a theorem due to the Horodeckis (7^ which can 
be stated as follows: any function B satisfying the condi- 
tions l)-3) below is an upper bound for the entanglement 
of distillation. 

1. Weak monotonicity: B(p) > B(A(p)) where A 
is any trace-preserving superoperator realizable by 
means of LOCC operations. 

2. Partial subadditivity: B{p n ) < nB(p). 

3. Continuity for isotropic states p p given by equation 
(55). Suppose that we have a sequence of isotropic 
states p p such that Tr{p p P + ) — > 1, if d — * oo. Then 
we require 



1 



lim 

d^oo l0g 2 d 



B( Pp 



1 



(72) 



We have already shown that LE n -i(p), for n > 1, sat- 
isfics conditions (1) and (2). From equation (54) 

LE n:1 (p p ) = log 2 (dp + i-^), V n > 1 (73) 

By evaluating this expression now for large d, we easily 
obtain that condition (3) is satisfied. □ 

It is also possible to state a proposition like (8) to 
the bounds on the distillable entanglement. E 1 will in 
this case provide an upper bound to Ep of all bipartite 
partitions. 



XII. LOWER BOUNDS FOR THE 
ENTANGLEMENT OF FORMATION 

One of the most celebrated entanglement measures is 
the entanglement of formation 



E F (p) 



mm piE E (\ipi}) 7 



(74) 



where Ee is the entropy of entanglement. Although this 
measure has a very meaningful physical interpretation 
and good properties, its calculation has been done only 
for a very class of states [78|. We show in this section 
that any entanglement witness can be used to provide 
lower bounds to the entanglement of formation. 

Let p — l'í) (í'l be a pure bipartite state with the fol- 
lowing Schimdt dccomposition: 

d 

\V)=J2 C MÏÏ> cr > c 2 > ... > c d . (75) 
i=i 

An analytic expression for the random robustness R r and 
the generalized robustness Rq of a pure state given by 
equation (78) is 0] 



r g {p)= \J2 c A - 1 ' 



(76) 
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R r (p) = cic 2 . 



(77) 



We start with two bounds for the entropy of entan- 
glement, i.e., the von Neumann entropy of the reduccd 
density matrix of a pure state \i Jj). In the case of two 
qubits, Wootters has shown that [TíJ 



H 



'l + yj\-\c{c( 



where H(x) — —x\og(x) — (1 — x) log(l — x). That is a 
particular case of the more general inequality 



H 



1 + Vi - 4c? 



2„2 



<-^ Cí 2 i g( C ?), J2 c i = 1 - 

i i 

(78) 



Anothcr similar inequality is 
Í>1 - 1 



log(d) - 1 



< 



5> 2 log(c 2 ). (79) 



Equations (72-73) can be proved maximizing the L.H.S. 
minus the R.H.S. and noting that the maximum is null 
in both cases. 

Choosing {pi, to be an optimal ensemble in equa- 
tion (71), we have 



Ef(p) = J2pM\^)) > Y,P' H ( 

i i \ 



> H 



1 + ^l-ARl{p) 



where we have used the convexity of R r and f(x) — 



H 



^ l+%/l-4x 2 ^ 



Similarly, we find 
log(d) - 1 



E F (p) > 



Rg(p). 



(80) 



The bound derived from R r is suitable for slightly entan- 
gled states, where the Schimdt coefhcients of the optimal 
\tjji) decay fast enough, making the truncation in the sec- 
ond Schimdt coefhcient a good approximation. 

As an first example, we consider the Horodecki 3x3 
states |80|. These states exhibit bound entanglement, 
since they have positive partial transposition. They are 
given by 



p{a) 



a 
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(81) 




FIG. 7: (Coloronline) Lower bound for the Horodecki states 
using R r . 



This family of states is interesting to test the first bound 
since their entanglement of formation was numerically 
calculated by Audenaert et al j8j| and was found to be 
very low. Figure (7) shows the bound provided by R r for 
the states p' = ep(a) + (1 - e)(I/D). 

For our second example we consider one of the unique 
states for which a analytic formula for Ep is know. It was 
shown by Terhal and Vollbrecht |Ï3 that for the isotropic 
states 82 1 



PF 



1-F 



7-P H 



FP 



E F ( PF )= d -^^l(F-l) + lo g ( d ), Fe 
From section (VII) we find 



4(d-l) 



d 2 



1 



E F ( PF ) > (log(d) - 1) [ F - - 



(82) 
(83) 



We see that, in this case, for sufficiently large d, the 
difference of the bound and the actual value of E F is 
always less than F. 

Note that every entanglement witness W, after being 
normalized such that either Tr(W) = 1 or W < I holds, 
can be used to deliver lower bounds to the entanglement 
of formation. 



XIII. ENTANGLEMENT, THERMODYNAMICS 
AND LATTICE SYSTEMS 

The study of entanglement properties of many-body 
systems, manly condensed matter, has received much 
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attention recently [Ü, [Ü, E3, Ü3, HE EE H3, [Ü, IÜ, 
l90l l9Ïl l92| . Several important models have been an- 
alyzed and connections with thermodynamic variables, 
such as internal energy and magnetization, have been 
raised [H ü . The negativity and concur- 
rence have been the most used measures, partly due their 
easy computation, but also because they made possible 
the derivation of some interesting simple thermodynamics 
like equations. This can be understood from the view of 
the witnessed entanglement. Every quantie derived from 
(7) not only defincs a measure of the degree of entangle- 
ment, but also gives a Hermitian operator, which vary for 
each state, whose expectation value quantifies the entan- 
glement of the state in question. It is exactly the possi- 
bility of measure experimentally the amount of entangle- 
ment, which is a feature shared by all common thermo- 
dynamics variables, that makes quantities expressed by 
(7) useful to the study of entanglement thermodynami- 
cal properties. In this section we present, as an example, 
the study of entanglement in the XXX Heisenberg model 
with and without a magnètic field using E^-i — Rq. The 
corresponding Hamiltonian is given by 



where the magnetization and the internal energy are 
given, respectively, by M = Yli (°í) anc ^ ^ = (-^0- 

We now proceed analysing the relation between en- 
tanglement and the magnètic susceptibility (x) m ther- 
mal states of i?2- According to Brukner, Vedral and 
Zeilinger 0, under temperatures close to zero and 
at zero external magnètic field, \ = (g 2 ^l 2 b / kT)[N I + 



2Ng 2 l· ij 3 2NR G 



3kT 



3kT 



(88) 



Remarkably, we see that the susceptibility is given by a 
term which resembles the classical Curie law more a term 
which takes into account the entanglement presented in 
the state. The equation above can be seen as a quanti- 
tative version of the experimental result of Ghose et al 
|92j | , who shown that at very low temperatures the mag- 
nètic susceptibility of certain materials is affected by the 
existence of entanglement. 



H 



XXX 



N 

»=i 



(84) 



i=l 



We first consider B = 0, in which case both Hamiltonians 
have SU(2) symmetry. According to section (VII), we 
can restrict the EWs in (7) to the ones that also have 
SU(2) symmetry. Then, from a Standard result from 
representation theory [g^, |m| , we find that all EWs with 
this symmetry can be written as 



(85) 



where Vi are unitary permutation operators. From an- 
alytic and numerical studies for the XXX Heisenberg 
model of odd N in the fundamental state and in the 
thermodynamical limit |95j |. we find that all other cor- 
relators are very small compared to the first neighbor 
correlators. We, thus, use the following ansatz for the 
optimal entanglement witness for the thermal states, at 
very low temperatures, of Hamiltonians (79) and (80) 



W =\ NI 



N 

E 

i=l 



/2N, 



where the factor 2N in the denominator comes from W < 
I. Note that this is the EW introduced by Toth et al [U 
Assuming that \B\ << | J|, and using the continuity of 
OEWs, we find that for the XXX Heisenberg model, at 
temperatures sufficiently close to zero, 



Rg 



U - BM 1 
2NJ 2' 



(87) 



XIV. CONCLUSION 



Summarizing, we have presented a new perspective to 
the quantification of entanglement based on witness op- 
erators. Several important EMs were shown to fit into 
this scenario and a new infinite family of EMs was in- 
troduced. The usefulness of the witnessed entanglement 
was illustrated by the study of diverse features of entan- 
glemnt, including super-selection rules constraints and 
efficiency of quantum information protocols. Finally, we 
have shown some interesting preliminary results in the 
study of thermodynamical properties of entanglement in 
macroscopic systems. 

We believe the results presented in this paper are only 
preliminary. The quantification on entanglement with 
EWs might be a very fruitful approach to development 
of the theory of entanglement, specially in the new appli- 
cations of entanglement, such as in identifying quantum 
phase transitions and improving the approximation of 
mean field theories 1971. 
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